High Frequency dynamics in metallic glasses 
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Using Inelastic X-ray Scattering we studied the collective dynamics of the glassy alloy Ni33Zre7 
in the first pseudo Brillouin zone, an energy-momentum region still unexplored in metallic glasses. 
We determine key properties such as the momentum transfer dependence of the sound velocity and 
of the acoustic damping, discussing the results in the general context of recently proposed pictures 
for acoustic dynamics in glasses. Specifically, we demonstrate the existence in this strong glass of 
well defined (in the Ioffe Regel sense) acoustic-like excitations well above the Boson Peak energy. 

PACS numbers: 61.43.Fs; 63.50.+x; 61.10.Eq 



The physics of disordered materials has attracted con- 
siderable interest during recent years [lj. This is mainly 
due to new theoretical concepts that can handle systems 
lacking long-range order, to large scale computer simula- 
tions based on realistic interactions and to the simultane- 
ous development of new experimental techniques. Disor- 
dered materials show several kinds of universal behaviors 
that differ from those of the ordered state. In particular, 
great interest has been devoted to their high frequency 
(THz) dynamics in respect to the corresponding crys- 
talline counterpart. In this context, Inelastic Neutron 
Scattering (INS) investigations of the collective dynam- 
ics in simple metallic glasses, which may be regarded as 
model systems, have played a leading role @> EI 0,01 Be- 
cause of their comparably simple atomic structure and 
their nearly isotropic interaction, metallic glasses (and 
glasses with Lennard- Jones (L-J) interaction), were the 
first tackled in extended computer simulations, and one- 
to-one comparisons with experiments were often under- 
taken. Simple systems like Mg7 Zn 3 o belong however to 
the minority of metallic glasses and interest has therefore 
moved soon to more abundant and characteristic repre- 
sentatives like the Zr-based transition metal-transition 
metal (TM-TM) or the TM-metalloid glasses, which are 
characterized by a strong topological and chemical short 
range order [|J. 

An important role in this context plays the TM-TM 
glass NiZr, which can be produced by melt spinning in 
a large range of compositions. Of these, the Zr-rich con- 
centrations, for which also crystalline NiZr^ exists, have 
played a prominent role in the investigation of the dy- 
namics of metallic glasses 0, 0- Concerning collective 
dynamics, in particular, the dispersion of modes around 
Q p , (the Q position of the main peak of the static struc- 
ture factor), i.e. in the region of maximal diffuse umk- 
lapp scattering §, was measured first for glassy Ni 2 4Zr7 6 
U and then in Ni67Zr33 [ljj. I n these pioneering exper- 
iments it was possible to prove the existence of Pseudo 
Brillouin Zones with clear Zone Boundaries near the max- 
ima of the static structure factor, these latter acting 



like "smeared out lattice points" 3]. However, impor- 
tant questions concerning the properties of the collec- 
tive excitations remained unanswered on the experimen- 
tal side: the existence of a low Q linear dispersion, the 
wave- vector and energy transfer dependence of the veloc- 
ity and the damping of the excitations. In this context, 
of importance here, is the recently proposed existence 
of a crossover frequency, which should mark the transi- 
tion between different dynamical regimes 0, ■ This 
crossover frequency could be identified in different ways. 
A first possibility (Qir) is by the Ioffe- Regel criterion, 
i.e. by the condition Sl/n = T, ft being the excitation 
frequency and T the sound attenuation. Furthermore, 
a crossover (f2 co ) could be associated to the frequency 
where the acoustic damping T (i.e. the broadening of the 
measured excitations) crosses from a T tx w 4 behavior 
(lo < fl co ) to a weaker power law, T cx u a with a ~ 2 
(lu > £! co ). Finally, the crossover frequency could be 
identified with the frequency of the Boson Peak (BP), 
i.e. the maximum in the excess of vibrational density 
of states g(u>) with respect to the Debye behavior of the 
corresponding crystal, represented as ■ Recently, the 
coincidence of these three different definition of crossover 
frequencies has been proposed to hold, at least for the 
case of strong glasses jlH . 

In this letter we report the determination of the col- 
lective dynamics in a metallic glass (M33Z1-67), focused 
to the region below Q p , where the excitations are still 
very well defined, aiming to answer some of the ques- 
tions on the nature of these high frequency modes, unan- 
swered up to now by INS investigations. Using Inelas- 
tic -X-ray Scattering (IXS), we exploited the lack of any 
kinematic restrictions in the accessible Q-E region which, 
combined with the good resolution, allowed us to investi- 
gate previously inaccessible dynamics in metallic glasses. 
A well denned acoustic branch has been observed up to 
Q ~ Q p , and the Ioffe- Regel criterium indicates that H,m 
is close to il(Q p ), i.e. well above the Boson Peak(BP) fre- 
quency. The behavior of the sound attenuation can be 
rationalized in terms of a relaxation process related to 
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the structural disorder, similarl y to what was found in 
simple monatomic systems |l3l Il4j. No evidence for a 
Rayleigh regime (T oc w 4 ) could be observed in the ex- 
plored Q range, thus concluding that httco, if existing, is 
well below 5 meV. Finally, the fragility value (m = 26), 
as deduced by the temperature dependence of the non 
ergodicity parameter, agrees well with independent vis- 
cosity determinations [15| , thus corroborating the recent 
proposed correlation between the low temperature vibra- 
tional regime and the high temperature diffusive dynam- 
ics in glass formers |l6|. and extending it to a class of 
systems not considered so far. The whole scenario re- 
ported here does not indicate the existence in this system 
of a single crossover frequency marking the BP position 
(TlUbp — 3 meV), the Ioffe-Regel criterium (hflm ~ 
8 meV) and the end of Rayleigh scattering of phonons 
regime (which, if it exists, has to be confined at energies 
Mlco < 5 meV). 

The experiment was performed at room temperature 
(T=297 K) at the IXS beam-line ID16 of ESRF H3, E3| 
at fixed scattering angle (and therefore fixed Q-values) 
in a region between 1.5 and 20 nm _1 . The Q-resolution 
(FWHM) was set to 5Q w 0.35 rnn" 1 and improved to 
5Q w 0.1 nm -1 at the lowest Q's. Using the (9 9 9) reflec- 
tion for the Si monochromator and crystal analyzers the 
overall energy resolution (FWHM) was 5%ui = 3.0 meV. 
A five-analyzers bench, operating in horizontal scatter- 
ing geometry, allowed the simultaneous collection of spec- 
tra at five different values of Q. Each energy scan from 
—50 < fiLO < 50 meV, where Ticu is the energy transfer 
(Eq — E), with Eq and E being the energy of the incident 
(17794 eV) and the scattered X-ray, took approximately 
300 minutes, and several scans have been summed up 
to improve the statistical accuracy. The Ni33Zr67 sam- 
ple was freshly prepared by melt spinning techniques, its 
structure was checked by X-ray diffraction before (and 
during) the measurement. 

In FigQ we report all the measured spectra in abso- 
lute counts, after proper normalization for the incident 
monitor. The elastic line is an effect of the non ergod- 
icity typical of the glassy state (scattering from static 
disorder) and possibly to incoherent scattering from con- 
centration fluctuations. Most reasonably the measured 
inelastic signal can be ascribed to a single acoustic prop- 
agating mode. The low energy {E w 8 meV) peak ob- 
served around Q ~ 11 nm _1 in the energy loss side is 
due to the elastic scattering from the kapton window 
located in the incoming beam before the sample. This 
peak, which is due to the sample environment geome- 
try and is only present in the range 9 < Q < 13 nm _1 , 
is also present in the empty cell measurements and to- 
tally disappears when the empty cell is subtracted from 
the data. The fitting procedure is performed by adding 
the empty cell with the proper transmission normaliza- 
tion factor to the model line-shape used to represent 
the data. Following generalized hydrodynamics, the sim- 
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FIG. 1: IXS spectra of Ni33Zr67 at the indicated fixed Q 
values (open dots). Also reported are the best-fit line-shapes 
(continuous lines, see text). The low energy structures at 
E ~ 8 meV observed around Q k 11 nm _1 are spurious 
scattering from the kapton windows of the sample container. 



plest model function to get the basic features of collec- 
tive glassy dynamics is the Damped Harmonic Oscillator 
(DHO), based on the assumption of an instantaneously 
decaying memory function for vibrational dynamics, typ- 
ical of Markovian processes, plus an clastic contribution 
due to the arrested diffusive dynamics and to incoherent 
scattering [19|. We have checked that fitting the data 
with more refined model (Debye-like memory functions), 
while improving the quality of the overall fit, does not 
produce significant variation in the peak width and po- 
sition. Therefore, for simplicity, we use the DHO model 
function: 



S(Q,u) 
S(Q) 



A(Q)6(w) 



l-A(Q) 



V 2 (Q)T(Q) 



7T (lj 2 -Q?(Q)) 2 +lj 2 T 2 (Q) 
(1) 

This function needs to be adapted to satisfy the de- 
tailed balance condition and to be convoluted with the 
instrumental energy resolution. The characteristic fre- 
quency of the acoustic mode, Q, corresponds to the max- 
imum of the longitudinal current spectra Cl{Q,u) = 
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(uj/Q) 2 S(Q,uj). Its Q dependency defines the disper- 
sion relation of the mode. The parameter T(Q) is related 
to the sound attenuation, while A(Q) is the intensity of 
the elastic scattering relative to the total intensity of the 
spectrum at the Q under consideration. 
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FIG. 2: Dispersion curve HQ(Q) of the NiZr alloy as obtained 
from the IXS and INS experiments. The momentum transfers 
have been scaled to Q p to account for the different structural 
properties due to the different compositions (Q p N * 33Zr67 ^ — 
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= 29 nm ). Inset: Q-dependence 



of the sound velocity v = D./Q as measured by IXS. 

Based on the form of the dispersion reported in Fig. [21 
we ascribe the corresponding excitations to longitudinal 
acoustic modes. In spite of the fact that we were able to 
measure at rather low Q-values, where the longitudinal 
acoustic and optic branch should be well separated and 
the excitations still well defined, we did not observe lon- 
gitudinal optic modes within the energy range covered 
in our experiment. This is most likely due to the fact 
that the intensity of the optic-like modes is expected to 
vanish in the Q — ► limit. Moreover, the dynamic struc- 
ture factor measured by the IXS is heavily dominated by 
the density-density correlations, while the concentration- 
concentration contribution is negligible: approximating 
the form factors of Ni and Zr by their Q — > val- 
ues, i.e. their atomic numbers, one has IixsiQi^) oc 
2.25S pp (Q,u) + 0.25S cc (Q,u) - 1.5S pc (Q,u). In the 
same figure, we report the energy of the high frequency 
mode previously measured with INS in Ni 6 7Zr 33 . Due 
to the different concentrations and cross-sections of the 
two species in the neutron and X-ray investigations, lead- 
ing to different weighting factors in the total dynamic 
structure factor measured in each of the two experiments, 
some difference in the measured excitation is observed. 
Nevertheless normalizing the momentum transfer scale 
to Q p , one obtains a good overlap. 

The sound attenuation, reported in Fig. 3a, shows a 



power law dependence, with a change of slope around 
Q m 7 nm -1 [E « 4 meV) being proportional to Q 2 
below this value (full line) and to a lower Q dependence 
above (dashed line). In the same figure we also report 
the peak position Sl(Q) divided by n (small full dots and 
thin full line). The intersection of this line with T(Q) de- 
fine the Ioffe-Regel limit, which is found here at Qir=\4 
±2 nm -1 or ±1 meV. The Boson peak energy, 

on the contrary, is Mlsp—i meV 0,H3|- This observa- 
tion clearly indicates that in the present glass there exist 
well defined (in the Ioffe-Regel sense) acoustic excitations 
well above the Boson peak energy. Furthermore, no tran- 
sition towards the Rayleigh Q 4 dependence is observed 
in the investigated Q range, thus posing an upper limit 
to the crossover energy, fiflco < 4.5meV. We can there- 
fore conclude that neither the Boson peak energy, nor the 
crossover energy, coincide with the Ioffe-Regel energy. 

The Q 2 dependence is a well documented feature of the 
high frequency attenuation in glasses, while the break- 
down of this dependence at intermediate Q (around Q—7 
nm -1 here) is not well understood up to now. In some 
systems, the breakdown of the quadratic dependence is 
simply an effect of the modulation of the structure of 
S(Q), which can be wiped out reporting the attenua- 
tion versus the excitation energy. In such cases, it seems 
more si gnifi cant to describe the sound attenuation as a 
f2 2 law [21|. This does not seem to be the case here, 
where the S(Q) is almost flat up to 10 nm -1 , while the 
change of slope is around 7 nm -1 , where the dispersion 
leaves its linear slope. Indeed, in Fig. Ob one can ob- 
serve a deviation from the r oc SI 2 law above fl ps 20 
meV (Q ~ 8 nm -1 ). The observed transition might be 
possibly related to a dynamical effect, as reported in the 
molecular dynamics simulation of a simple metallic glass 
(open triangles in fig Et) , scaled by an arbitrary factor 
for presentation reasons) [14) . 

In Fig. the elastic/total intensity ratio A(Q) = ^ sL - 
is reported. Generally speaking, the elastic signal comes 
from the relaxation processes active over a timescale 
longer than the inverse of the phonon frequency. In 
systems, where the dominant elastic contribution comes 
from the structural relaxation, A(Q) is a measure of the 
so called non ergodicity factor, f(Q). It has been re- 
cently proposed [l6| that the low temperature behavior 
of the non ergodicity factor in a T/T g plot is closely 
correlated to the so called Angel plot, i.e. to the tem- 
perature behavior of viscosity on approaching the glass 
transition in a T g /T scale. More specifically, in terms of 

the fragility index m = \iuiT^T g d(T S )r) anc ^ °^ ^ ne P a ~ 
rameter a = limT^o dfi^W^y one nas m ~ 135a. At the 
same time, the temperature dependence of the non ergod- 
icity parameter in the T — > limit is well described by 
the functional form f(Q, T — > 0) = 1+ J r / T ■ The glass 
transition temperature of Ni33Zr67 is about T g — 652 K, 
and using our determination of A(Q — > 0,T) w 0.92 at 
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FIG. 3: (a): Q dependency of the damping of the observed 
excitation (full triangles). A transition from the Q 2 to a lin- 
ear dependence of the damping around Q = 7 nm _1 can be 
observed. Values of 0,/ir are also reported to identify the 
Ioffe-Regel crossover (full circles), together with the Boson 
Peak energy. Open triangles show the sound attenuation cal- 
culated by means of molecular dynamics in a simple metal, 
which shows a qualitatively similar behavior |l4|. (b): Energy 
dependence of the damping: the breakdown of the quadratic 
dependence is shifted to higher energies but not completely 
removed, (c): reduced VDOS ^ for Ni33Zr67 as measured 
by INS [23 and calculated by MD simulations 0. 
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FIG. 4: Q-dependence of the elastic/total scattered intensity. 
Inset: static structure factor, note the prepeak at Q ~ 14 
nm -1 ), where A(Q) changes its slope. 

possible transition towards a Q 4 dependence. Finally, 
we determine the Ioffe-Regel limit, that is reached close 
to the top of the acoustic branch, at Ejp=7 ±1 meV, 
thus well above the Boson peak energy, Ebp= 3 meV. In 
the present strong glass, there exist well defined (in the 
Ioffc Regel sense) acoustic-like excitations well above the 
Boson Peak energy. These findings are not in accordance 
with recent mod els p ro p osed for the acoustic properties 
in strong glasses [Tl|,ll2| that predict Ejr^Ebp- 
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room temperature one has a — 0.19, which would lead 
to a fragility of about m = 26, indicating Ni33Zr67 to be 
a strong glass. Though this value of fragility has to be 
taken as a lower estimate, due to possible elastic contri- 
bution from concentration fluctuations, it is in very good 
agreement with independent high temperature viscosity 
determination, which would give m = 24 15] . 

Summing up, we have studied the high frequency 
dynamics of the strong metallic glass M33Z1-67 in a 
range of momentum transfer below the first maximum 
of the static structure factor. Evidence of a longitudi- 
nal acoustic-like branch has been reported, which is the 
extension into the sensible low Q region of the high fre- 
quency excitations observed a few years ago by means 
of INS. The analysis of the non-ergodicity factor allows 
us to classify -in agreement with previous determination 
[po| - this glass as "strong" (m—26). At low Q, the sound 
attenuation exhibits a Q 2 behavior, with no hint of a 
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